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ABSTRACT 

Supernovae type la (SNIa) are used as “standard candles” for cosmological distance 
scales. To fit their light curve shape - absolute luminosity relation, one needs to assume 
an intrinsic color and a likelihood of host galaxy extinction or a convolution of these, 
a color distribution prior. The host galaxy extinction prior is typically assumed to be 
an exponential drop-off for the current supernova programs {P{Ay) oc We 

explore the validity of this prior using the distribution of extinction values inferred 
when two galaxies accidentally overlap (an occulting galaxy pair). We correct the 
supernova luminosity distances from the SDSS-III Supernova projects (SDSS-SN) by 
matching the host galaxies to one of three templates from occulting galaxy pairs based 
on the host galaxy mass and the Ay-bias - prior-scale (tq) relation from Jha et al. 
(2007). 

We find that introducing an Ay prior that depends on host mass results in lowered 
luminosity distances for the SDSS-SN on average but it does not reduce the scatter 
in individual measurements. This points, in our view, to the need for many more 
occulting galaxy templates to match to SNIa host galaxies to rule out this possible 
source of scatter in the SNIa distance measurements. We match occulting galaxy 
templates based on both mass and projected radius and we find that one should 
match by stellar mass first with radius as a secondary consideration. 

We discuss the caveats of the current approach: the lack of enough radial coverage, 
the small sample of priors (occulting pairs with HST data), the effect of gravitation¬ 
ally interacting as well as occulting pairs, and whether an exponential distribution is 
appropriate. Our aim is to convince the reader that a library of occulting galaxy pairs 
observed with HST will provide sufficient priors to improve (optical) SNIa measure¬ 
ments to the next required accuracy in Cosmology. 

Key words: 


1 INTRODUCTION 

Supernova Type la (SNIa) are a prime Cosmological tool 
(Riess et al. 1998; Perlmutter et al. 1999). Their relation 
between light curve and intrinsic luminosity make them ex¬ 
cellent standard candles. Part of the exquisite precision re¬ 
quired is a good understanding of the expected dust extinc¬ 
tion and known reddening-dimming relation, the Extinction 
Law (Cardelli et al. 1989; Calzetti et al. 1994). With the 
Dark Energy detection result in hand, the supernova com- 
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munity now looks for the next step in accuracy (e.g., the 
3% solution Riess et al. 2011). A much better model for 
the dust extinction in the Supernova’s host galaxy will be 
critical. Host galaxy dust is already identified by the Dark 
Energy Task Force as a principal unknown (Albrecht et al. 
2006) and this remains the case, together with the photo¬ 
metric uncertainty across surveys (e.g., Conley et al. 2011; 
Scolnic et al. 2014; Betoule et al. 2014). 

These issues are being addressed by the current gen¬ 
eration of supernovae searches, e.g., the Supernova Legacy 
Survey (SNLS, Astier et al. 2006), the Equation of State: 
SupErNovae trace Cosmic Expansion (ESSENCE, Wood- 
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Vasey et al. 2007), and the Sloan Digital Sky Survey Su- 
perNova survey (SDSS-SN, Kessler et al. 2009; Sako et al. 
2014). The supernovae community has been systematically 
exploring what could affect the observed relation between 
peak luminosity and light curve width (Phillips 1993) for 
SNela. 

A prime candidate is the correlation with properties of 
the galaxy hosting the SNIa (Gallagher et al. 2005; Sullivan 
et al. 2006; Mannucci et al. 2006; Pan et al. 2014; Kistler 
et al. 2013). Two populations of SNIa appear associated 
with two populations of stars: blue, star forming galaxies 
host higher rates of fast declining supernova and red, passive 
galaxies host predominantly more slowly-declining SNe fa 
(Hamuy et al. 1996; Howell 2001; van den Bergh et al. 2005; 
Mannucci et al. 2006; Schawinski 2009; Lampeitl et al. 2010; 
Wang et al. 2013). The difference in light curve characteris¬ 
tics may be completely attributable to the different progen¬ 
itor stellar populations. 

However, simultaneously, the issue of host galaxy ex¬ 
tinction has come to the fore, either the applicability of the 
appropriate extinction law (Riess et al. 1996; Phillips et al. 
1999; Altavilla et al. 2004; Reindl et al. 2005; Jha et al. 
2007; Conley et al. 2007), the validity of the prior of ex¬ 
tinction values (Jha et al. 2007; Wood-Vasey et al. 2007; 
Gupta et al. 2011), or the possibility that the geometry of 
the dusty ISM in (and hence the extinction distribution) 
may change with galaxy evolution (Holwerda 2008; Holw¬ 
erda et al. 2015). Host galaxies appear to be indistinguish¬ 
able from normal field galaxies (Childress et al. 2013). After 
photometric calibration, this host galaxy extinction issue is 
likely to return as a dominant issue standing in the way of 
the community’s final goal of 1% precision for SNI (e.g., 
Riess et al. 2011; Kelly et al. 2014). 

There is now an opportunity to probe the dependen¬ 
cies of the extinction prior through two projects. First, the 
SDSS-SN project has generated a uniform sample of SNIa 
and their host galaxy properties of exquisite quality. Sec¬ 
ondly, occulting galaxy pairs promise to map extinction in 
spiral disks using the Hubble Space Telescope (HST) imag¬ 
ing, enabling a way to generate an extinction prior for SN la 
fit. The first of these extinction maps are now available. 

In this paper, we explore the effects of applying the 
extinction priors from occulting galaxy pairs to the SDSS- 
SN sample and gauge the effects. We opt for the SDSS-SN 
sample among all the supernova samples available now as 
it has uniform properties of not only the supernovae, as 
measured with both light curve packages, but also on the 
host galaxies. The latter is important if one wants to cor¬ 
rect for the effects of host galaxy properties. The SDSS-SN 
project uses two light curve fitting packages, MLCS2 k 2 and 
SALT2. We use the former’s output as it uses an explicit 
prior for host galaxy extinction, making a direct comparison 
with the occulting galaxy technique practical, even though 
SALT and its derivative BaSALT are now often the more 
preferred package (e.g., Guy et al. 2010; Betoule et al. 2014; 
Rest et al. 2014). We should note here that the Ay value 
used in MLCS2 k 2 is not the exact same quantity as mea¬ 
sured from galaxy transparency measurements, i.e., occult¬ 
ing galaxy measurements. MLCS2 k2’s Ay parameter is the 
difference in supernova color from an assumed intrinsic one, 
translated to an extinction via an assumed reddening law 
(in the case of the SDSS-SN Kessler et al. 2009, a grayer 



Figure 1. The 2MASX occulting galaxy pair at z=0.06. The ex¬ 
tinction in the overlap region (black aperture) can be estimated 
from the complementary apertures; the foreground spiral (F, the 
green aperture) and the background galaxy (B, red aperture). 


one than Milky Way). The reddening law in this instant 
is an implicit input setting, one which has been tweaked to 
minimize Gosmological fit residuals (e.g. Hicken et al. 2009). 

We assume in this paper that the shape of the Ay fit 
parameter’s prior and the shape of the distribution of Ay 
values observed in occulting galaxies is to first order the 
same. 

We explore the extinction effects by matching SN la 
host galaxies to those occulting galaxies with HST data 
through their stellar mass. In a previous paper on SN la (?), 
we explored the relation between the extinction distribution 
and inclination and found that an inclination correction us¬ 
ing a simple cos(i) is sufficient to make the Ay distribution 
found by the SDSS-SN identical in high and low-inclination 
subsamples. We therefore adopt this inclination correction 
in this paper as well. We will discuss the caveats and limits 
of the current approach and sketch out how to proceed from 
here in the generation of improved extinction priors, tailored 
to the host galaxy of the observed supernova. 


2 OCCULTING GALAXIES 

One can accurately measure the extinction of light by in¬ 
terstellar dust using a known light source. The most accu¬ 
rate extra-galactic method is to use a partially occulting 
galaxy pair (Figure 1). Assuming both galaxies are symmet¬ 
ric, one can estimate flux contributions to the overlap region 
from the complementary regions at the same radius in each 
galaxy. The reliance on differential photometry enables an 
extinction measurement {Ay) in a single-color image with¬ 
out relying on any assumed Extinction Law. With multiple 
filters, the effective Extinction Law itself can be measured 
(Keel & White 2001a; Holwerda et al. 2009; Keel et al. 2014). 


2.1 The Occulting Galaxies Method 

Figure 1 shows an example pair, serendipitously imaged 
with HST for the ANGST survey (Dalcanton et al. 2009), 
near NGG 253 (2MASX) To determine the extinction in the 
overlap region (black aperture), one estimates the rela¬ 
tive flux contributions by both galaxies from the comple¬ 
mentary apertures in the foreground (green aperture) and 
background galaxy (red aperture in Figure 1), respectively, 
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assuming rotational symmetry (i.e. the image be can ro¬ 
tated and be self-similar). Thus, we estimate the flux con¬ 
tributions to the overlap region from the same image. In 
the case of ground-based images, each of these apertures 
would yield single Ay value, averaged over a large portion 
of the foreground disk. In the case of an HST image, it re¬ 
sults in a highly accurate map of the dust extinction in the 
overlap region with hundreds of independent lines-of-sight. 
The map is constructed pixel-by-pixel; for each pixel in the 
overlap aperture, the corresponding pixel in the background 
galaxy aperture (red) is found i.e., the pixel at the same 
distance from the background galaxy center but in the op¬ 
posite direction, as well as the corresponding pixel in the 
foreground galaxy (same distance with respect to the fore¬ 
ground galaxy center). We estimate the optical depth (r) or 
extinction {Ay) from the flux in the overlap {F -\- Be~^), 
a mix of flux from background and foreground galaxies, the 
flux in the corresponding background pixel (B) and fore¬ 
ground pixel {F, same colors as Figure 1): 


Thus, a map of dust extinction and a distribution func¬ 
tion {P{Ay)) such as Figure 2 can be constructed from an 
HST image. For each type of foreground galaxy, inclination 
and radial distance, an appropriate P{Ay) can thus be ob¬ 
served. 


2.2 Occulting Pairs 

The accuracy in individual pairs is limited by the symmetry 
of both galaxies and image quality. The uncertainty from 
the assumption of galaxy symmetry can mitigated by using 
elliptical galaxies as the background galaxy and by combin¬ 
ing the results for several pairs. In the case of HST images 
(e.g.. Figure 1), a single pair contains enough independent 
lines-of-sight through the foreground galaxy to construct a 
P{Ay) such as Figure 2. 

The occulting galaxies technique however has increased 
steadily in accuracy and usefulness, owing in a large part 
to the increasing sample sizes. Estimating dust extinction 
and mass from differential photometry in occulting pairs of 
galaxies was first proposed by White & Keel (1992). Their 
technique was then applied to all known pairs using ground- 
based optical images (Andredakis & van der Kruit 1992; 
Berlind et al. 1997; Domingue et al. 1999a; White et al. 
2000) and spectroscopy (Domingue et al. 2000), and later 
space-based HST images (Keel & White 2001a,b; Elmegreen 
et al. 2001; Holwerda et al. 2009, 2011). These initial results 
were limited by sample sizes. More recently, more pairs were 
found in the SDSS spectroscopic catalog (86 pairs in Hol¬ 
werda et al. 2007), and in an ongoing effort in the Galaxy- 
ZOO project (Lintott et al. 2008). The current count stands 
at 1993 pairs (Keel et al. 2013) providing opportunities for 
follow-up with IFU observations (Holwerda et al. 2013; Hol¬ 
werda & Keel 2013) and GALEX (Keel et al. 2014). 

Highly complete spectroscopic surveys, such the Galaxy 
and Mass Assembly (GAMA Driver et al. 2009, 2011; Baldry 
et al. 2010) are adding many more overlapping pairs through 
the identification of blended spectra or close on-sky pairs 
(Holwerda et al. in prep.). Observational results include 



Extinction (Ay) 

Figure 2. The shape of the Ay distribution of the three occult¬ 
ing galaxies with enough overlap to generate a distribution. These 
distributions have not yet been corrected for their foreground 
disk’s inclination. Dashed lines are exponential fits through the 
distributions. The exponential drop-off can then be used as a 
SNIa prior. The derived drop-offs are listed in Figure 5, with 
their tq values listed in Table 1. 

mean extinction profiles (White et al. 2000; Domingue et al. 
2000; Holwerda et al. 2007), an indication that the dust 
may be fractal (Keel & White 2001a) and that the color- 
extinction relation is gray in the ground-based observations 
Holwerda et al. (2007). The Galactic Extinction Law re¬ 
turns as soon as the physical sampling of the overlap re¬ 
gion resolves the molecular clouds in the foreground disk 
(< lOOpc, Keel & White 2001a,b; Elmegreen et al. 2001; 
Holwerda et al. 2009). 

The extinction maps of foreground disks made with 
HST are especially useful as a probability map of extinction 
towards a single object of interest. To date, we have obtained 
HST imaging of three overlapping pairs with a large enough 
overlap region to generate a reasonable distribution. The 
normalized distributions of these three are shown in Figure 
2. Radial coverage is very similar to the radii at which SNIa 
are found (Table 1): between 0 and 6 effective radii (con¬ 
verted using the prescriptions in Graham & Driver 2005). 

Our previous results on occulting galaxy pairs are con¬ 
sistent with an exponential decline with radius of the average 
disk attenuation, not dissimilar to the exponential profile of 
the stellar light (e.g., Domingue et al. 1999b; Holwerda 2005; 
Holwerda et al. 2005a,b). We stress here that we deal here 
now in the exponential distribution of values in a certain 
range of radii but this is not the same as the radial expo¬ 
nential decline. For example, the lower values {Ay ~ 0.1) are 
typically, but not exclusively, found at higher radius while 
the higher values are more typical for spiral arms. 

2.3 UGC 3995 

This pair of galaxies has been long known to be overlapping 
and interacting (e.g., Marziani et al. 1999) and we have an¬ 
alyzed the extinction properties of both CALIFA DRl IFU 
data and archival WFPC2 F606W imaging in Holwerda & 
Keel (2013). Radial coverage is from 5-15”, covering the full 


Ay — 1.086 X r = —1.086 x In 


{F + Be- 


B 
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Table 1. The properties of the foreground occulting galaxy in 
three pairs with sufficient HST information to form an extinction 
map. Stellar mass in solar masses and radii in effective (half-light) 
radii. 


Pair 

Type 

Stellar Mass 
M* 

Radial Coverage 

(^e) 

(-fo) 

UGC 3995 

Sa 

1.7 X 10i° 

0-2 

0.22 

2MASX 

Sd 

3.1 X 10® 

2-5 

0.85 

AM 0500 

Sbc 

2.6 X 10i° 

~ 3.7 

0.67 



Figure 3. The distribution of SNla host galaxy mass from Sako 
et al. (2014). They present both FSPS (Conroy & Gunn 2010) 
and PEGASE (Bruzual A. 2009) stellar masses. The masses of 
2MASX (thick dashed), UGC 3995 (dashed line), and AM500 
(dotted line) are marked as well. 


15”60 Petrosian radius (SDSS-r) of UGC3995B, the fore¬ 
ground galaxy. We obtained the foreground galaxy’s mass 
from the SDSS fluxes of the unocculting side and the pre¬ 
scription from Zibetti et al. (2009). 


2.4 2MASX 

This pair’s serendipitous discovery was originally reported 
in Holwerda et al. (2009). The distance and stellar mass was 
estimated in Holwerda et al. (2013) and the radial coverage 
of the overlap region is between 2 and 5 effective (or half- 
light) radii {Re). 


2.5 AM 500-620 

This was a pair identified in White et al. (2000) as an over¬ 
lapping galaxy pair and the HST data were originally re¬ 
ported in Keel & White (2001a). The radial coverage is 
mostly in the outer disk (~ i? 25 , the 25 mag/arcsec^ as 
defined by de Vaucouleurs et al. 1991). Stellar mass and 
Hubble type are from Keel & White (2001a). 


3 SUPERNOVAE TYPE lA SAMPLE 

Type la Supernovae (SNIa) are the touchstone of Cosmo¬ 
logical distance measurements. The photometric properties 
of their light curves are very stable but appear to depend 
slightly on host galaxy properties (e.g., star-formation rate 
or stellar mass). With this mind, the third incarnation of the 
Sloan Digital Sky Survey (SDSS-HI) included a supernova 
search with spectroscopic follow-up. This SDSS-Supernova 
search (SDSS-SN) ^ has yielded a wealth of galaxy and su¬ 
pernova properties. 

SDSS-SN is described in Frieman et al. (2008) and in 
more detail in Sako et al. (2008) and Sako et al. (2014). 
Observing for three months of the year from 2005 to 2007, 
SDSS-SN identified hundreds spectroscopically confirmed 
SNe la in the redshift range 0.05 < z < 0.35. 

Sako et al. (2014) present the hnal data-product of this 
tremendous observational effort. They include an assessment 
of SN type based on the light curve (PSNID output) and 
light curve fits using the two most commonly used packages, 
SALT2 (Guy et al. 2007) and MLCS2 k 2 (Jha et al. 2007). 
Of these two packages, the SALT2 light curve fitter is now 
the most popular. Yet we focus on the MLCS2 k 2 package 
because it explicitly starts with a host galaxy extinction 
prior, making it easier to directly compare its assumptions 
to our observations in occulting pairs. 

However, we should note that both SALT2 and 
MLCS2 k 2 essentially measure the reddening of the super¬ 
nova, not actual attenuation but the convolution of intrinsic 
supernova color distribution, scatter by the surrounding in¬ 
terstellar matter and dust attenuation in the host galaxy. 
The shape of the MLCS2 k2 Av prior used for the light 
curve fits was: P(Av) ~ with tq = 0.4. Because it 

has host galaxy extinction as an explicit prior and a model 
relation between Ay distribution width (tq) and Ay bias, 
we use the MLCS2 k 2 values for our further analysis here. 
The MLCS2 k 2 version used for the SDSS-SN analysis as¬ 
sumed a Ry = 2.8 reddening relation to convert reddening 
to extinction Ay (Kessler et al. 2009), indicating it is not 
completely host attenuation alone. However, it is not as un¬ 
physical as the Ry = 1.7 needed elsewhere with MLCS2 k 2 
(Hicken et al. 2009). For our purposes here, we will explic¬ 
itly assume that the shape of the prior for MLCS2 k2’s ht 
parameter Ay can be taken to be similar to the distribution 
of actual attenuation values as found in occulting galaxies. 

Galaxy properties are those available from the SDSS- 
DR9 database and stellar mass and star-formation are mod¬ 
eled with two different packages, FSPS (Gonroy & Gunn 
2010) and PEGASE (Bruzual A. 2009). For the purpose 
of this paper, the stellar mass is of interest and we choose 
the FSPS value for for the further analysis. Figure 3 shows 
the distribution of host mass galaxies with SN la. The mean 
mass is ~ lO^^M©, similar to AM500. We focus only on those 
objects that have reasonable chance of being bona-hde SN la 
(S/N>5 and P(SNIa) > 50%) to ensure the conclusions for 
the improved prior are based on these only (1698 SN in the 
total sample). 

In addition to the values thoughtfully provided by the 
SDSS-SN project in Sako et al. (2014), we retrieve the axis 
ratio, position angle, and Petrosian radius from the SDSS 

^ http://www.sdss.org/supernova/aboutsupernova.html 
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Figure 4. The distribution of the projected separation of SN la 
from their host galaxy in effective radii (Re) based on the values 
reported in Sako et al. (2014). For comparison, Sako et al. (2014) 
plot the distance-to-host or "directional light radius” (do^ij), 
which is not deprojected into the plane of the host disk. 


server. Using these values, we compute the galactocentric 
radial distance from the center of the host galaxy in Pet¬ 
rosian radii. This is a different value than that of the sepa¬ 
ration presented in Sako et al. (2014), the ’’directional light 
radius” {(Idlr), which is not deprojected into the plane of 
the host disk. If we assume disks (n=l Sersic index), the 
Petrosian radius corresponds to approximately two effective 
radii {Rp — 2.15 x Re, see Graham & Driver 2005). Figure 
4 shows the radial coverage for the SNIa for which SDSS 
Petrosian radii were available. SNIa are typically found by 
SDSS-SN between 0 and 5 effective radii (< 2.5 Petrosian 
radii), which correspond approximately to the radii covered 
by the occulting galaxy pairs (Table 1). 

4 ANALYSIS 

4.1 Av Distribution and Bias 

The relation between the exponential drop-off and bias in¬ 
troduced in the Av value was explored by Jha et al. (2007) 
using the Riello & Patat (2005) simulations. Figure 5 shows 
their relation (Figure 22 in the appendix of Jha et al. 2007). 
We mark the ro values inferred from the occulting galaxy 
pairs, corrected for inclination (to = T/cos{i). To parame¬ 
terize the Av bias, we fit the relation in Figure 5 with: 

AA = -0.5 X -b 0.02 x ro, (2) 

where ro is the exponential drop-off of the host galaxy prior 
under consideration. 

We note first that none of our three occulting galaxy 
templates match the drop-off typically used in MLCS2K2 
(to = 0.4) which was used for the SDSS-SN analysis (Sako 
et al. 2008, 2014). Therefore, assuming the SDSS-SN galax¬ 
ies adhere to one of these three templates based on their host 
mass will inevitably introduce a bias in the distribution. Our 
goal is to explore in which direction and how severe it would 
be. 



Figure 5. The relation between prior exponential scale and the 
resulting bias in the Av from Jha et al. (2007). The thick gray 
line is the MLCS2K2 prior. The values of the three occulting pairs 
are also marked (not corrected for their inclination). 



Redshift (z) 


Figure 6. The relation between redshift and MLCS2K2 luminos¬ 
ity distance (top panel) and the difference between the ACDM 
and the SNIa luminosity distances (bottom panel). The relation 
for different equations of state (ju = 1/ — 1) are shown for illus¬ 
tration only. 


4.2 New Luminosity Distances for the SDSS-SN 

Figure 6 shows the relation between the luminosity distance 
and redshift for the SDSS-SN Snia. They closely follow the 
expected ACDM relation with scatter around the mean in 
the residual (bottom panel in Figure 6). To extract the ex¬ 
act nature of Dark Energy, one approach is to explore this 
residual. However, inconsistent or heterogeneous photome¬ 
try or host galaxy extinction could introduce biases in this 
residual. 

We assign each host galaxy a new exponential template 
{toco) from one of the three occulting galaxy (OCG) distri¬ 
butions (fit as an exponential drop-off), based on the nearest 
in stellar mass, and correct for disk inclination: 
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Figure 7. The distribution of host galaxy r^ost values as calcu¬ 
lated from equation 3. There is still a peak at tq = 0.35 as this is 
the value assigned when no good host estimate was available. The 
concentration at tq = 0. is due to poorly constrained inclination. 


Thost = 'rocG(M) X cos{ih ost ) 


( 3 ) 


where ihost is the inclination of the host galaxy as inferred 
from the minor and major axes ratio: 


COs{ihost) — 


{h/aY - (&/a)i 


1 - ib/a) 


2 

min 


( 4 ) 


originally from Hubble (1926), where a is the major axis, b 
the minor axis and we assume bjamin = 0.1. We found in 
Holwerda et al. {submitted) that the distribution of SDSS-SN 
Av values becomes the same for high- and low-inclination 
galaxies, once the values are corrected for disk inclination 
using this simple cos{i) correction (i.e., the screen approxi¬ 
mation). 

Figure 7 shows the distribution of resulting Thost for the 
SDSS-SN SNIa sample. The peak around Thost = 0.4 corre¬ 
sponds to all the host galaxies for which no mass estimate 
was available in the SDSS-SN catalog. 

Starting from the new Thost exponential values for each 
host galaxy (equation 3), we can now calculate the Ay bias 
for each host, based on its stellar mass and inclination (equa¬ 
tion 2) and apply this to the luminosity distance. Figure 8 
shows the relative difference with ACDM distances as a func¬ 
tion of redshift and the mean values (and dispersion) before 
our Av bias correction and after (Table 2). 

There are two things to note in the corrected values. 
First is that the mean luminosity distance is slightly lower, 
i.e., there is a little more extinction in the SN la host galaxies 
if we use individual priors, but there is no less scatter in the 
values. 

The change in the mean is interesting for cosmological 
study but at present only points to the direction one can 
expect the luminosity distances to change. The priors may 
be slightly more tailored but they are hardly appropriate 
for each host galaxy (see below for the occulting galaxy pair 
prior caveats). This lack of an appropriate prior is reflected 
-in our opinion- in the lack of change in dispersion. If an 
appropriate correction to the prior shape and hence Ay bias 


Figure 8. The relative difference between ACDM predicted dis¬ 
tance and the SN la luminosity distances from SDSS-SN. Indi¬ 
vidual SNIa (grey points) and the mean and rms (red points) 
and once corrected for Ay bias (green points) using the occulting 
galaxy templates (Figure 2), matched through host mass. Alter¬ 
natively, one can match OCG templates by radius (blue points). 
Two Cosmologies are shown as well to illustrate: w=-0.9 and -1, 
neither a fit to the data. The different Cosmologies are shown to 
illustrate how the choice of dust prior can have an almost similar 
effect on the inferred Cosmology. 


Table 2. The average and scatter for A^/D^ in four bins for 
the original MLCS2K2 analysis, one corrected using the occulting 
galaxy templates matched with host mass or radius. 


z 

Az,/Dl 

Mass-Matched 

Radius Matched 

0.10 

0.32 (0.60) 

0.29 (0.61) 

-0.05 (0.79) 

0.20 

0.13 (0.24) 

0.11 (0.24) 

0.06 (0.34) 

0.30 

0.10 (0.25) 

0.08 (0.25) 

0.04 (0.34) 

0.40 

0.17 (0.32) 

0.15 (0.32) 

0.13 (0.37) 


had been applied, luminosity distances would adhere closer 
to the Hubble flow and any remaining scatter can be at¬ 
tributed to individual motion of the galaxies or photometric 
errors. 

Alternatively, one could match the few occulting galaxy 
templates by their projected radius, not the foreground 
galaxy’s mass. These are the blue averages points in Fig¬ 
ure 8. The scatter around the new mean residuals is worse 
than in the case of the mass-matched templates. We can 
speculate at this point that to match templates, one should 
select by host mass first and galactocentric radius only as a 
secondary effect. Once again, we are restricted by the num¬ 
bers of occulting pair templates we have. 


5 CAVEATS 

While the radial coverage and host mass range appears rea¬ 
sonably approximated by our three template host galaxies, 
there are a number of issues outstanding with just using 
three galaxies as templates. We will focus on the shortcom¬ 
ings of the occulting galaxy pairs presented as priors. 
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Table 3. The exponentially modified Gaussian fits to the Ay 
distributions observed in occulting pairs. 


pair 

m 

a 

A 

2MASX 

0.013 

0.015 

6.65 

UGG3559 

0.092 

0.092 

2.10 

AM500 

0.013 

0.011 

2.29 

AM1316 

1.41 

0.50 

70.60 


5.1 MLCS2K2 Ay parameter is all Dust 
Attenuation 

The central assumption in this work has been that the at¬ 
tenuation seen in occulting galaxies is the same as the Ay 
parameter in MLCS2K2. However, because the Ay is de¬ 
rived assuming an attenuation law, it cannot be ruled out 
that some measure of the reddening is not from dust atten¬ 
uation but scatter in, or SN intrinsic color scatter. 

5.2 Radial Coverage 

While Table 1 and Figure 4 seem to suggest that at least the 
SN la occur mostly at the same radii for which we have some 
extinction distribution from occulting pairs, the coverage is 
far from complete. For example, for low-mass galaxies, there 
is only information between 1-2 effective radii. 

5.3 Small Sample 

It is obvious to point out that it is impossible to infer a re¬ 
lation between stellar mass and an extinction distribution 
with only three galaxy templates. We plan to obtain more 
templates from future HST observations (Holwerda 2014) 
but for now we can only assign the exponential prior to the 
closest stellar mass to the host mass. Alternatively, we fit¬ 
ted a linear relation between tq and mass for the occulting 
galaxy pairs but this resulted in no improvement of the vari¬ 
ance of the SN la distances. 

5.4 Interactions 

Occulting galaxies are often interacting and the redshift dif¬ 
ference for two of our templates (UGC 3995 and 2MASX) 
are such that an ongoing interaction cannot be ruled out. It 
is very likely in the case of UGC 3995. The exact effect of 
gravitational interaction on the distribution of ISM is un¬ 
known but likely to be severe. Is there more dust now at 
higher radii? Do the shocks and tides move more dust into 
denser clouds in advance of the burst of star-formation asso¬ 
ciated with gas-rich interactions? The way around this issue 
again is to observe more and bona-fide occulting pairs, well 
separated in redshift. A selection from a spectroscopic sur¬ 
vey such as the Galaxy And Mass Assembly Survey (GAMA 
Driver et al. 2009, 2011; Baldry et al. 2010) may well prove 
to be ideal. 

5.5 Is the Ay prior really an exponential? 

One strong assumption so far has been that the distribution 
of extinction values is an exponential drop-off. However, Fig- 



Extinction (Ay) 

Figure 9. The shape of the Ay distribution of the three oc¬ 
culting galaxies with enough overlap to generate a distribution. 
An exponentially modified Gaussian distribution was with using 
maximum likelihood to each distribution. The relevant parame¬ 
ters are listed in Table 3. 

ure 2 shows that none of the three distributions is fully de¬ 
scribed as an exponential. For example, in the more massive 
galaxies, the peak of the distribution is not close to Ay = 0 
at all. We noted this in Holwerda & Keel (2013) for both 
the HST observations and the IFU ones. In the case of UGC 
3995, one could contribute this to the ongoing interaction. 
Again the only way forward it to obtain more templates for 
the same host mass and infer the mean and variance of the 
distribution function. 

Our assumption -and the MLCS2K2 authors’- is that 
the shape of the distribution is an exponential. We now 
strongly suspect the true shape of the extinction distribu¬ 
tions may be a log-normal one, similar to what Dalcanton et 
al. (in preparation) find for the extinction values in front of 
stars in M31 based on the PHAT project (Dalcanton et al. 
2012) or an exponentially modified Gaussian. To illustrate, 
Figure 9 shows fits to the Ay distribution using the expo¬ 
nentially modified Gaussian fit to each distribution using 
maximum likelihood. 

This strongly suggests that MLCS2K2 will now have to 
be redone using these templates instead of an exponential 
one. 


6 APPLICATION TO SALT2 

SALT2 handles the SNIa color and dust extinction in a 
combined parameter, and it may not be obvious how our 
approach has any easy application to this package. To illus¬ 
trate, Figure 10 shows the relation between the observed Ay 
and the SALT2 color parameter, c (similar to Figure 9 in 
Sako et al. 2008). There is a linear relation between the two 
that follows: 

c = 0.046 * Av -b 0.012 (5) 

Thus, a change in the prior for Ay in the case of MLCS2 k 2 
would translate a change in the prior for SALT2’s assump¬ 
tion but in color-space. To illustrate in figure 11, we plot the 
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Figure 10. The relation between extinction as measured by the 
MLCS2 k 2 package and the color parameter from the SALT2 
package for the SDSS-SN supernovae. There is a linear relation be¬ 
tween the two (fit with maximum likelihood, including the SALT2 
color errors. 



Figure 11. The distribution of SALT2 color parameter c for 
all SN la (black line) and after converting the Ay values from 
MLCS2 k 2 using equation 5 (green histogram). Since the fit is 
reasonable, there is no big difference between the two. After cor¬ 
rection for inclination alone, the Ay values would translate into 
the blue histogram. 


distribution of c values found in SDSS-III for SN la, as well 
as the values of Ay from MLCS2 k 2 translated using the 
above relation. The two distributions are -unsurprisingly- 
very similar, however, if we then translate the Ay value cor¬ 
rected for the inclination of the disk (inferred from the axis 
ratio), the distribution of c values changes significantly. In 
our view, this change in distribution after inclination correc¬ 
tion is the original distribution one could adopt as the prior 
distribution of the parameter c for SALT2 applications, sim¬ 
ilar to the one to be adopted for MLCS2 k 2 (since both of 
them are essentially color priors). 


7 DISCUSSION 

The SNIa community is well aware of the effects of host 
galaxy extinction may have on the results and have ac¬ 
tively been searching for both “dust-free” supernovae (i.e., 
hosted by early-type galaxies, see e.g., Suzuki et al. 2012) 
and pushed light-curve observations into the near-infrared 
to mitigate extinction issues. 

This paper’s aim is an attempt to map out how the 
available observational information, and specifically the host 
galaxy information, can be further used to optimize the 
SNIa distance measurement. As a first exercise, we ap¬ 
plied the three available templates from occulting galax¬ 
ies, naively matched to the SNIa host through the stellar 
mass or galactocenric radius. Introducing new host extinc¬ 
tion templates certainly has an effect bnt since the scatter 
in the distances is not markedly reduced, using just these 
templates would just skew the Cosmological result, not im¬ 
prove it. Despite the fact that the Ay values were attained 
using a more Milky-Way type relation (Kessler et al. 2009), 
we cannot rule out that none of the current scatter comes 
from host attenuation and all of it comes from intrinsic color 
scatter in the SN la. But assuming much of the scatter is due 
to host attenuation, when matching host templates through 
projected radius rather than host mass, the scatter increases 
more than through mass-matching. This points, in our opin¬ 
ion, in which order one should match the occulting galaxy 
templates; first throngh mass and then through galactocen- 
tric radius. 

Another key issue for SN la cosmological measurements 
will be if the extinction prior will change with redshift. Given 
that between z=0 and z=l, specific star-formation rates in¬ 
creased by a factor 3-10 (e.g., Noeske et al. 2007), one can ex¬ 
pect a mnch more fractured dusty ISM (thanks to increased 
turbulence, as we speculated in Holwerda 2008) or a much 
denser one (to sustain the star-formation). Some initial re¬ 
sults on dust lanes in L* edge-ons suggest that the global 
morphology in disk galaxies remains the same: 80% of these 
have a dust lane mid-plane (Holwerda et al. 2012b). 

Alternatively, if the Ay distribution is significantly dif¬ 
ferent for massive versus lower-mass disk galaxies, any evo¬ 
lution in the mix of SNIa host galaxies will result in an 
effective change in the appropriate Ay prior. Dust in is 
strongly suspected to be distributed differently with respect 
to the stars in lower-mass galaxies. Initial Spitzer observa¬ 
tions showed very extended dust morphology in a few indi¬ 
vidual cases of lower-mass galaxies (Hinz et al. 2007, 2009), 
similar to the occulting pair in Figure 1. Recently, it has be¬ 
come clear that in the stacked SED of lower-mass galaxies 
there is virtually no attenuation in the optical and nltra- 
violet yet Herschel emission clearly show large dust masses 
at very low temperatures (Popescu et al. 2002; Gross! et al. 
2010; Bourne et al. 2012). The combined evidence points to a 
dependence of star/dust geometry on galaxy mass. Whether 
this means dust is relatively more extended in radius, height 
above the plane (e.g., Dalcanton et al. 2004; Holwerda et al. 
2012a) or dumpiness, remains to be seen. 

Fortunately, global extinction measures can be checked 
using high-redshift occulting galaxy pairs and both effects 
can then be acconnted for, be it evolution in the mix of 
host galaxies or in the ISM of the host galaxies themselves. 
There is, at presence, no observed evolution in the color 
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distribution of observed SNIa (Conley et al. 2011; Betoule 
et al. 2014). However, if there is significant evolution in the 
P{Av) distribution with redshifts -either due to ISM evolu¬ 
tion in the host galaxies or a significant change in the host 
galaxy population mix-, it may change the relative lumi¬ 
nosity distances as a function of redshift, the driver or dark 
energy cosmology. 

We hope to have convinced the reader that the use of 
occulting galaxy pairs to obtain extinction distribution tem¬ 
plates is a good way forward to help reduce the remaining 
uncertainties in SNIa distance measurements. The eventual 
goal is to cover enough range in stellar mass, radial cover¬ 
age and inclinations to map the mean and variance of the 
typical extinction distribution in spiral galaxy disks. This 
will have many astrophysical applications but the use as a 
prior for Cosmological SN la measurements was the principal 
driver for us to obtain HST observations of occulting pairs: 
GO-13695, STarlight Absorption Reduction through a Sur¬ 
vey of Multiple Occulting Galaxies (STARSMOG) P.I. B.W. 
Holwerda, a 150 orbit SNAPshot program (Holwerda 2014), 
which is projected to complete a sample of 98 occulting pairs 
by 2017. 


ACKNOWLEDGEMENTS 

The authors are very grateful for the insightful comments 
of the anonymous referee and to the SDSS-SN team who 
provided the data. KJM is supported by an Australian Re¬ 
search Council Discovery Early Career Researcher Award. 
This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propul¬ 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Admin¬ 
istration. This research has made use of NASA’s Astro¬ 
physics Data System. This research made use of Astropy, 
a community-developed core Python package for Astron¬ 
omy (Astropy Collaboration et al. 2013). This research made 
use of matplotlib, a Python library for publication quality 
graphics (Hunter 2007). PyRAF is a product of the Space 
Telescope Science Institute, which is operated by AURA for 
NASA. This research made use of SciPy (Jones et al. 2001). 


REFERENCES 

Albrecht A., Bernstein G., Cahn R., Freedman W. L., He¬ 
witt J., Hu W., Huth J., Kamionkowski M., Kolb E. W., 
Knox L., Mather J. C., Staggs S., Suntzeff N. B., 2006, 
astro-ph /0609591 

Altavilla G., Fiorentino G., Marconi M., Musella I., Cap- 
pellaro E., Barbon R., Benetti S., Pastorello A., Riello M., 
Turatto M., Zampieri L., 2004, MNRAS, 349, 1344 

Andredakis Y. C., van der Kruit P. C., 1992, A&A, 265, 
396 

Astier P., Guy J., Regnault N., Pain R., Aubourg E., Balam 
D., Basa S., Carlberg R. G., Fabbro S., Fouchez D., Hook 
I. M., Howell D. A., Lafoux H., Neill J. D., Palanque- 
Delabrouille N., Perrett K., Pritchet C. J., Rich J., Sulli¬ 
van M., Taillet R., Aldering G., Antilogus P., Arsenijevic 
V., Balland C., Baumont S., Bronder J., Courtois H., Ellis 
R. S., Filiol M., Gongalves A. G., Goobar A., Guide D., 


Hardin D., Lusset V., Lidman G., McMahon R., Mouchet 

M. , Mourao A., Perlmutter S., Ripoche P., Tao G., Walton 

N. , 2006, A&A, 447, 31 

Astropy Collaboration, Robitaille T. P., Tollerud E. J., 
Greenfield P., Droettboom M., Bray E., Aldcroft T., Davis 

M. , Ginsburg A., Price-Whelan A. M., Kerzendorf W. E., 
Gonley A., Grighton N., Barbary K., Muna D., Fergu¬ 
son H., Grollier F., Parikh M. M., Nair P. H., Unther 
H. M., Deil C., Woillez J., Gonseil S., Kramer R., Turner 
J. E. H., Singer L., Fox R., Weaver B. A., Zabalza V., 
Edwards Z. L, Azalee Bostroem K., Burke D. J., Casey 

A. R., Crawford S. M., Dencheva N., Ely J., Jenness T., 
Labrie K., Lim P. L., Pierfederici F., Pontzen A., Ptak A., 
Refsdal B., Servillat M., Streicher O., 2013, A&A, 558, 
A33 

Baldry 1. K., Robotham A. S. G., Hill D. T., Driver S. R, 
Liske J., Norberg P., Bamford S. P., Hopkins A. M., Love- 
day J., Peacock J. A., Gameron E., Groom S. M., Cross 

N. J. G., Doyle I. F., Dye S., Frenk G. S., Jones D. H., 
van Kampen E., Kelvin L. S., Nichol R. C., Parkinson 
H. R., Popescu G. G., Prescott M., Sharp R. G., Suther¬ 
land W. J., Thomas D., Tuffs R. J., 2010, MNRAS, 404, 
86 

Berlind A. A., Quillen A. C., Pogge R. W., Sellgren K., 
1997, AJ, 114, 107 

Betoule M., Kessler R., Guy J., Mosher J., Hardin D., 
Biswas R., Astier P., El-Hage P., Konig M., Kuhlmann 

S., Marriner J., Pain R., Regnault N., Balland C., Bassett 

B. A., Brown P. J., Gampbell H., Carlberg R. G., Gellier- 
Holzem F., Cinabro D., Conley A., D’Andrea C. B., De- 
Poy D. L., Doi M., Ellis R. S., Fabbro S., Filippenko 
A. V., Foley R. J., Frieman J. A., Fouchez D., Galbany 
L., Goobar A., Gupta R. R., Hill G. J., Hlozek R., Hogan 
G. J., Hook I. M., Howell D. A., Jha S. W., Le Guillou 

L. , Leloudas G., Lidman C., Marshall J. L., Moller A., 
Mourao A. M., Neveu J., Nichol R., Olmstead M. D., 
Palanque-Delabrouille N., Perlmutter S., Prieto J. L., 
Pritchet G. J., Richmond M., Riess A. G., Ruhlmann- 
Kleider V., Sako M., Schahmaneche K., Schneider D. P., 
Smith M., Sollerman J., Sullivan M., Walton N. A., 
Wheeler C. J., 2014, A&A, 568, A22 

Bourne N., Maddox S. J., Dunne L., Auld R., Baes M., 
Baldry I. K., Bonheld D. G., Gooray A., Groom S. M., 
Dariush A., de Zotti G., Driver S. P., Dye S., Eales 

S., Gomez H. L., Gonzalez-Nuevo J., Hopkins A. M., 
Ibar E., Jarvis M. J., Lapi A., Madore B., Michalowski 

M. J., Pohlen M., Popescu G. G., Rigby E. E., Seib¬ 
ert M., Smith D. J. B., Tuffs R. J., van der Werf P., 
Brough S., Buttiglione S., Cava A., Clements D. L., Con- 
selice C. J., Fritz J., Hopwood R., Ivison R. J., Jones 
D. H., Kelvin L. S., Liske J., Loveday J., Norberg P., 
Robotham A. S. G., Rodighiero G., Temi P., 2012, ArXiv 
e-prints/1201.1916 

Bruzual A. G., 2009, A&A, 500, 521 

Galzetti D., Kinney A. L., Storchi-Bergmann T., 1994, 
ApJ, 429, 582 

Gardelli J. A., Clayton G. C., Mathis J. S., 1989, ApJ, 345, 
245 

Ghildress M., Aldering G., Antilogus P., Aragon C., Bailey 

S., Baltay G., Bongard S., Buton C., Canto A., Cellier- 
Holzem F., Chotard N., Copin Y., Fakhouri H. K., Gan- 
gler E., Guy J., Hsiao E. Y., Kerschhaggl M., Kim A. G., 


© 2010 RAS, MNRAS 000, 1-?? 



10 B. W. Holwerda 


Kowalski M., Loken S., Nugent P., Paech K., Pain R., 
Pecontal E., Pereira R., Perlmutter S., Rabinowitz D., 
Rigault M., Runge K., Scalzo R., Smadja G., Tao C., 
Thomas R. C., Weaver B. A., Wu C., 2013, ApJ, 770, 
107 

Conley A., Carlberg R. G., Guy J., Howell D. A., Jha S., 
Riess A. G., Sullivan M., 2007, ApJ, 664, L13 

Conley A., Guy J., Sullivan M., Regnault N., Astier P., Bal- 
land C., Basa S., Carlberg R. G., Fouchez D., Hardin D., 
Hook I. M., Howell D. A., Pain R., Palanque-Delabrouille 
N., Perrett K. M., Pritchet C. J., Rich J., Ruhlmann- 
Kleider V., Balam D., Baumont S., Ellis R. S., Fabbro S., 
Fakhouri H. K., Fourmanoit N., Gonzalez-Gaitan S., Gra¬ 
ham M. L., Hudson M. J., Hsiao E., Kronborg T., Lidman 

C. , Mourao A. M., Neill J. D., Perlmutter S., Ripoche P., 
Suzuki N., Walker E. S., 2011, ApJS, 192, 1 

Conroy C., Gunn J. E., 2010, ApJ, 712, 833 

Dalcanton J. J., Williams B. F., Lang D., Lauer T. R., 
Kalirai J. S., Seth A. C., Dolphin A., Rosenfield P., Weisz 

D. R., Bell E. F., Bianchi L. C., Boyer M. L., Caldwell 
N., Dong H., Dorman C. E., Gilbert K. M., Girardi L., 
Gogarten S. M., Gordon K. D., Guhathakurta P., Hodge 
P. W., Holtzman J. A., Johnson L. C., Larsen S. S., Lewis 
A., Melbourne J. L., Olsen K. A. G., Rix H.-W., Rosema 

K. , Saha A., Sarajedini A., Skillman E. D., Stanek K. Z., 
2012, ApJS, 200, 18 

Dalcanton J. J., Williams B. F., Seth A. C., Dolphin A., 
Holtzman J., Rosema K., Skillman E. D., Cole A., Girardi 

L. , Gogarten S. M., Karachentsev 1. D., Olsen K., Weisz 
D., Ghristensen C., Freeman K., Gilbert K., Gallart C., 
Harris J., Hodge P., de Jong R. S., Karachentseva V., Ma¬ 
teo M., Stetson P. B., Tavarez M., Zaritsky D., Governato 
F., Quinn T., 2009, ApJS, 183, 67 

Dalcanton J. J., Yoachim P., Bernstein R. A., 2004, ApJ, 
608, 189 

de Vaucouleurs G., de Vaucouleurs A., Corwin H. G., Buta 

R. J., Paturel G., Fouque P., 1991, Third Reference Gat- 
alogue of Bright Galaxies. Volume 1-3, XII, 2069 pp. 7 
figs.. Springer-Verlag Berlin Heidelberg New York 

Domingue D. L., Keel W. C., Ryder S. D., White HI R. E., 
1999a, AJ, 118, 1542 

—, 1999b, AJ, 118, 1542 

Domingue D. L., Keel W. C., White HI R. E., 2000, ApJ, 
545, 171 

Driver S. R, Hill D. T., Kelvin L. S., Robotham A. S. G., 
Liske J., Norberg P., Baldry 1. K., Bamford S. P., Hop¬ 
kins A. M., Loveday J., Peacock J. A., Andrae E., Bland- 
Hawthorn J., Brough S., Brown M. J. L, Cameron E., 
Ching J. H. Y., Colless M., Conselice C. J., Croom S. M., 
Cross N. J. G., de Propris R., Dye S., Drinkwater M. J., 
Ellis S., Graham A. W., Grootes M. W., Gunawardhana 

M. , Jones D. H., van Kampen E., Maraston C., Nichol 
R. C., Parkinson H. R., Phillipps S., Pimbblet K., Popescu 
C. G., Prescott M., Roseboom 1. G., Sadler E. M., Sansom 
A. E., Sharp R. G., Smith D. J. B., Taylor E., Thomas D., 
Tuffs R. J., Wijesinghe D., Dunne L., Frenk G. S., Jarvis 
M. J., Madore B. F., Meyer M. J., Seibert M., Staveley- 
Smith L., Sutherland W. J., Warren S. J., 2011, MNRAS, 
413, 971 

Driver S. P., Norberg P., Baldry 1. K., Bamford S. P., Hop¬ 
kins A. M., Liske J., Loveday J., Peacock J. A., Hill D. T., 
Kelvin L. S., Robotham A. S. G., Cross N. J. G., Parkin¬ 


son H. R., Prescott M., Conselice C. J., Dunne L., Brough 

S. , Jones H., Sharp R. G., van Kampen E., Oliver S., Rose- 
boom 1. G., Bland-Hawthorn J., Croom S. M., Ellis S., 
Cameron E., Cole S., Frenk C. S., Couch W. J., Alister 
W. G., Proctor R., De Propris R., Doyle 1. F., Edmondson 

E. M., Nichol R. G., Thomas D., Eales S. A., Jarvis M. J., 
Kuijken K., Lahav O., Madore B. F., Seibert M., Meyer 
M. J., Staveley-Smith L., Phillipps S., Popescu C. C., San¬ 
som A. E., Sutherland W. J., Tuffs R. J., Warren S. J., 
2009, Astronomy and Geophysics, 50, 050000 

Elmegreen D. M., Kaufman M., Elmegreen B. G., Brinks 
E., Struck C., Klaric M., Thomasson M., 2001, AJ, 121, 
182 

Frieman J. A., Bassett B., Becker A., Choi C., Cinabro 
D., DeJongh F., Depoy D. L., Dilday B., Doi M., Gar- 
navich P. M., Hogan C. J., Holtzman J., Im M., Jha S., 
Kessler R., Konishi K., Lampeitl H., Marriner J., Mar¬ 
shall J. L., McGinnis D., Miknaitis G., Nichol R. C., Pri¬ 
eto J. L., Riess A. G., Richmond M. W., Romani R., 
Sako M., Schneider D. P., Smith M., Takanashi N., Tokita 
K., van der Hey den K., Yasuda N., Zheng C., Adelman- 
McGarthy J., Annis J., Asset R. J., Barentine J., Bender 

R. , Blandford R. D., Boroski W. N., Bremer M., Brewing- 
ton H., Collins C. A., Crotts A., Dembicky J., Eastman 
J., Edge A., Edmondson E., Elson E., Eyler M. E., Fil- 
ippenko A. V., Foley R. J., Frank S., Goobar A., Gueth 

T. , Gunn J. E., Harvanek M., Hopp U., Ihara Y., Ivezic 
Z., Kahn S., Kaplan J., Kent S., Ketzeback W., Kleinman 

S. J., Kollatschny W., Kron R. G., Krzesihski J., Lament! 
D., Leloudas G., Lin H., Long D. C., Lucey J., Lup- 
ton R. H., Malanushenko E., Malanushenko V., McMillan 
R. J., Mendez J., Morgan G. W., Morokuma T., Nitta A., 
Ostman L., Pan K., Rockosi C. M., Romer A. K., Ruiz- 
Lapuente P., Saurage G., Schlesinger K., Snedden S. A., 
Sollerman J., Stoughton C., Stritzinger M., Subba Rao 
M., Tucker D., Vaisanen P., Watson L. G., Watters S., 
Wheeler J. C., Yanny B., York D., 2008, AJ, 135, 338 

Gallagher J. S., Garnavich P. M., Berlind P., Ghallis P., 
Jha S., Kirshner R. P., 2005, ApJ, 634, 210 

Graham A. W., Driver S. P., 2005, Publications of the As¬ 
tronomical Society of Australia, 22, 118 

Gross! M., Hunt L. K., Madden S., Vlahakis G., Bomans 
D. J., Baes M., Bendo G. J., Bianchi S., Boselli A., 
Clemens M., Corbelli E., Cortese L., Dariush A., Davies 
J. L, de Looze L, di Serego Alighieri S., Fadda D., Fritz 
J., Garcia-Appadoo D. A., Gavazzi G., Giovanardi C., 
Hughes T. M., Jones A. P., Pierini D., Pohlen M., Saba- 
tini S., Smith M. W. L., Verstappen J., Xilouris E. M., 
Zibetti S., 2010, A&A, 518, L52-f 

Gupta R. R., D’Andrea C. B., Sako M., Conroy C., Smith 
M., Bassett B., Frieman J. A., Garnavich P. M., Jha S. W., 
Kessler R., Lampeitl H., Marriner J., Nichol R. C., Schnei¬ 
der D. P., 2011, ApJ, 740, 92 

Guy J., Astier P., Baumont S., Hardin D., Pain R., Reg¬ 
nault N., Basa S., Garlberg R. G., Conley A., Fabbro S., 
Fouchez D., Hook 1. M., Howell D. A., Perrett K., Pritchet 
C. J., Rich J., Sullivan M., Antilogus P., Aubourg E., 
Bazin G., Bronder J., Filiol M., Palanque-Delabrouille N., 
Ripoche P., Ruhlmann-Kleider V., 2007, A&A, 466, 11 

Guy J., Sullivan M., Conley A., Regnault N., Astier P., Bal- 
land C., Basa S., Carlberg R. G., Fouchez D., Hardin D., 
Hook 1. M., Howell D. A., Pain R., Palanque-Delabrouille 


© 2010 RAS, MNRAS 000, 1-?? 



Matching SN la Dust Extinction to a host 11 


N., Perrett K. M., Pritchet C. J., Rich J., Ruhlmann- 
Kleider V., Balam D., Baumont S., Ellis R. S., Fabbro S., 
Fakhouri H. K., Fourmanoit N., Gonzalez-Gaitan S., Gra¬ 
ham M. L., Hsiao E., Kronborg T., Lidman C., Mourao 
A. M., Perlmutter S., Ripoche P., Suzuki N., Walker E. S., 
2010, A&A, 523, A7 

Hamuy M., Phillips M. M., Suntzeff N. B., Schommer R. A., 
Maza J., Smith R. G., Lira P., Aviles R., 1996, AJ, 112, 
2438 

Hicken M., Wood-Vasey W. M., Blondin S., Ghallis R, Jha 

S. , Kelly P. L., Rest A., Kirshner R. P., 2009, ApJ, 700, 
1097 

Hinz J. L., Engelbracht C. W., Willmer C. N. A., Rieke 
G. H., Rieke M. J., Smith P. S., Blaylock M., Gordon 
K. D., 2009, in The Evolving ISM in the Milky Way and 
Nearby Galaxies 

Hinz J. L., Rieke M. J., Rieke G. H., Willmer G. N. A., 
Misselt K., Engelbracht G. W., Blaylock M., Pickering 

T. E., 2007, ApJ, 663, 895 

Holwerda B., 2014, STarlight Absorption Reduction 
through a Survey of Multiple Occulting Galaxies 
(STARSMOG). HST Proposal, SNAP-13695, 104 orbits 
Holwerda B. W., 2005, PhD thesis, Proefschrift, Rijksuni- 
versiteit Groningen, 2005 
—, 2008, MNRAS, 386, 475 

Holwerda B. W., Bianchi S., Boker T., Radburn-Smith D., 
de Jong R. S., Baes M., van der Kruit P. G., Xilouris M., 
Gordon K. D., Dalcanton J. J., 2012a, A&A, 541, L5 
Holwerda B. W., Boker T., Dalcanton J. J., Keel W. G., 
de Jong R. S., 2013, MNRAS, 433, 47 
Holwerda B. W., Dalcanton J. J., Keel W. G., Boeker T., 
de Jong R. S., 2011, MNRAS 
Holwerda B. W., Dalcanton J. J., Radburn-Smith D., de 
Jong R. S., Guhathakurta P., Koekemoer A., Allen R. J., 
Boker T., 2012b, ApJ, 753, 25 
Holwerda B. W., Gonzalez R. A., Allen R. J., van der Kruit 
P. C., 2005a, AJ, 129, 1381 
—, 2005b, AJ, 129, 1396 

Holwerda B. W., Keel W. C., 2013, A&A, 556, A42 
Holwerda B. W., Keel W. G., Bolton A., 2007, AJ, 134, 
2385 

Holwerda B. W., Keel W. G., Williams B., Dalcanton J. J., 
de Jong R. S., 2009, AJ, 137, 3000 
Holwerda B. W., Reynolds A., Smith M., Kraan-Korteweg 
R. G., 2015, MNRAS, 446, 3768 
Howell D. A., 2001, ApJ, 554, L193 
Hubble E. R, 1926, ApJ, 64, 321 

Hunter J. D., 2007, Computing In Science & Engineering, 
9, 90 

Jha S., Riess A. G., Kirshner R. P., 2007, ApJ, 659, 122 
Jones E., Oliphant T., Peterson P., Others, 2001, SciPy: 

Open source scientific tools for python 
Keel W. C., Manning A. M., Holwerda B. W., Lintott C. J., 
Schawinski K., 2014, AJ, 147, 44 
Keel W. C., Manning A. M., Holwerda B. W., Mezzoprete 
M., Lintott C. J., Schawinski K., Gay P., Masters K. L., 
2013, PASP, 125, 2 

Keel W. C., White HI R. E., 2001a, AJ, 121, 1442 
—, 2001b, AJ, 122, 1369 

Kelly P. L., Filippenko A. V., Burke D. L., Hicken M., 
Ganeshalingam M., Zheng W., 2014, ArXiv e-prints 
Kessler R., Becker A. C., Cinabro D., Vanderplas J., Frie- 


man J. A., Marriner J., Davis T. M., Dilday B., Holtzman 

J. , Jha S. W., Lampeitl H., Sako M., Smith M., Zheng C., 
Nichol R. C., Bassett B., Bender R., Depoy D. L., Doi 
M., Elson E., Filippenko A. V., Foley R. J., Garnavich 
P. M., Hopp U., Ihara Y., Ketzeback W., Kollatschny W., 
Konishi K., Marshall J. L., McMillan R. J., Miknaitis G., 
Morokuma T., Mortsell E., Pan K., Prieto J. L., Rich¬ 
mond M. W., Riess A. G., Romani R., Schneider D. P., 
Sollerman J., Takanashi N., Tokita K., van der Heyden 

K. , Wheeler J. C., Yasuda N., York D., 2009, ApJS, 185, 
32 

Kistler M. D., Stanek K. Z., Kochanek C. S., Prieto J. L., 
Thompson T. A., 2013, ApJ, 770, 88 
Lampeitl H., Smith M., Nichol R. C., Bassett B., Cinabro 
D., Dilday B., Foley R. J., Frieman J. A., Garnavich P. M., 
Goobar A., Im M., Jha S. W., Marriner J., Miquel R., 
Nordin J., Ostman L., Riess A. G., Sako M., Schneider 
D. P., Sollerman J., Stritzinger M., 2010, ApJ, 722, 566 
Lintott C. J., Schawinski K., Slosar A., Land K., Bamford 
S., Thomas D., Raddick M. J., Nichol R. C., Szalay A., 
Andreescu D., Murray P., Vandenberg J., 2008, MNRAS, 
389, 1179 

Mannucci F., Della Valle M., Panagia N., 2006, MNRAS, 
370, 773 

Marziani P., D’Onofrio M., Dultzin-Hacyan D., Sulentic 
J. W., 1999, AJ, 117, 2736 

Noeske K. G., Weiner B. J., Faber S. M., Papovich C., 
Koo D. C., Somerville R. S., Bundy K., Conselice C. J., 
Newman J. A., Schiminovich D., Le Floc’h E., Coil A. L., 
Rieke G. H., Lotz J. M., Primack J. R., Barmby P., Cooper 

M. C., Davis M., Ellis R. S., Fazio G. G., Guhathakurta 
P., Huang J., Kassin S. A., Martin D. C., Phillips A. C., 
Rich R. M., Small T. A., Willmer C. N. A., Wilson G., 
2007, ApJ, 660, L43 

Pan Y.-C., Sullivan M., Maguire K., Hook 1. M., Nugent 
P. E., Howell D. A., Arcavi L, Botyanszki J., Cenko S. B., 
DeRose J., Fakhouri H. K., Gal-Yam A., Hsiao E., Kulka- 
rni S. R., Laher R. R., Lidman C., Nordin J., Walker E. S., 
Xu D., 2014, MNRAS, 438, 1391 
Perlmutter S., Aldering G., Goldhaber G., Knop R. A., Nu¬ 
gent P., Castro P. G., Deustua S., Fabbro S., Goobar A., 
Groom D. E., Hook 1. M., Kim A. G., Kim M. Y., Lee 
J. C., Nunes N. J., Pain R., Pennypacker C. R., Quimby 
R., Lidman C., Ellis R. S., Irwin M., McMahon R. G., 
Ruiz-Lapuente P., Walton N., Schaefer B., Boyle B. J., 
Filippenko A. V., Matheson T., Fruchter A. S., Panagia 

N. , Newberg H. J. M., Couch W. J., The Supernova Cos¬ 
mology Project, 1999, ApJ, 517, 565 

Phillips M. M., 1993, ApJ, 413, L105 
Phillips M. M., Lira P., Suntzeff N. B., Schommer R. A., 
Hamuy M., Maza J., 1999, AJ, 118, 1766 
Popescu C. C., Tuffs R. J., Volk H. J., Pierini D., Madore 

B. F., 2002, ApJ, 567, 221 

Reindl B., Tammann G. A., Sandage A., Saha A., 2005, 
ApJ, 624, 532 

Rest A., Scolnic D., Foley R. J., Huber M. E., Chornock 
R., Narayan G., Tonry J. L., Berger E., Soderberg A. M., 
Stubbs C. W., Riess A., Kirshner R. P., Smartt S. J., 
Schlafly E., Rodney S., Botticella M. T., Brout D., Ghallis 
P., Czekala L, Drout M., Hudson M. J., Kotak R., Leibler 

C. , Lunnan R., Marion G. H., McCrum M., Milisavljevic 

D. , Pastorello A., Sanders N. E., Smith K., Stafford E., 


© 2010 RAS, MNRAS 000, 1-?? 



12 B. W. Holwerda 


Thilker D., Valenti S., Wood-Vasey W. M., Zheng Z., Bnr- 
gett W. S., Chambers K. C., Denneau L., Draper P. W., 
Flewelling H., Hodapp K. W., Kaiser N., Kudritzki R.- 
P., Magnier E. A., Metcalfe N., Price P. A., Sweeney W., 
Wainscoat R., Waters C., 2014, ApJ, 795, 44 

Riello M., Patat F., 2005, MNRAS, 362, 671 

Riess A. G., Filippenko A. V., Challis P., Clocchiatti A., 
Diercks A., Garnavich P. M., Gilliland R. L., Hogan C. J., 
Jha S., Kirshner R. P., Leibundgnt B., Phillips M. M., 
Reiss D., Schmidt B. P., Schommer R. A., Smith R. C., 
Spyromilio J., Stnbbs C., Suntzeff N. B., Tonry J., 1998, 
AJ, 116, 1009 

Riess A. G., Macri L., Casertano S., Lampeitl H., Fergnson 
H. G., Filippenko A. V., Jha S. W., Li W., Chornock R., 
2011, ApJ, 730, 119 

Riess A. G., Press W. H., Kirshner R. P., 1996, ApJ, 473, 
588 

Sako M., Bassett B., Becker A., Cinabro D., DeJongh F., 
Depoy D. L., Dilday B., Doi M., Frieman J. A., Garnavich 
P. M., Hogan C. J., Holtzman J., Jha S., Kessler R., Kon- 
ishi K., Lampeitl H., Marriner J., Miknaitis G., Nichol 
R. G., Prieto J. L., Riess A. G., Richmond M. W., Romani 

R. , Schneider D. P., Smith M., Subba Rao M., Takanashi 
N., Tokita K., van der Heyden K., Yasuda N., Zheng 
C., Barentine J., Brewington H., Ghoi C., Dembicky J., 
Harnavek M., Ihara Y., Im M., Ketzeback W., Klein- 
man S. J., Krzesihski J., Long D. C., Malanushenko E., 
Malanushenko V., McMillan R. J., Morokuma T., Nitta 
A., Pan K., Saurage G., Snedden S. A., 2008, AJ, 135, 348 

Sako M., Bassett B., Becker A. C., Brown P. J., Campbell 
H., Cane R., Cinabro D., D’Andrea C. B., Dawson K. S., 
DeJongh F., Depoy D. L., Dilday B., Doi M., Filippenko 
A. V., Fischer J. A., Foley R. J., Frieman J. A., Galbany 

L., Garnavich P. M., Goobar A., Gnpta R. R., Hill G. J., 
Hayden B. T., Hlozek R., Holtzman J. A., Hopp U., Jha 

S. W., Kessler R., Kollatschny W., Leloudas G., Marriner 
J., Marshall J. L., Miquel R., Morokuma T., Mosher J., 
Nichol R. C., Nordin J., Olmstead M. D., Ostman L., 
Prieto J. L., Richmond M., Romani R. W., Sollerman J., 
Stritzinger M., Schneider D. P., Smith M., Wheeler J. C., 
Yasuda N., Zheng C., 2014, ArXiv e-prints 

Schawinski K., 2009, MNRAS, 397, 717 

Scolnic D., Rest A., Riess A., Huber M. E., Foley R. J., 
Brout D., Chornock R., Narayan G., Tonry J. L., Berger 
E., Soderberg A. M., Stubbs G. W., Kirshner R. P., Rod¬ 
ney S., Smartt S. J., Schlafly E., Botticella M. T., Ghallis 
P., Czekala I., Drout M., Hudson M. J., Kotak R., Leibler 

C. , Lunnan R., Marion G. H., McGrum M., Milisavljevic 

D. , Pastorello A., Sanders N. E., Smith K., Stafford E., 
Thilker D., Valenti S., Wood-Vasey W. M., Zheng Z., Bur- 
gett W. S., Chambers K. C., Denneau L., Draper P. W., 
Flewelling H., Hodapp K. W., Kaiser N., Kudritzki R.- 
P., Magnier E. A., Metcalfe N., Price P. A., Sweeney W., 
Wainscoat R., Waters C., 2014, ApJ, 795, 45 

Sullivan M., Le Borgne D., Pritchet C. J., Hodsman A., 
Neill J. D., Howell D. A., Carlberg R. G., Astier P., 
Aubourg E., Balam D., Basa S., Gonley A., Fabbro 
S., Fouchez D., Guy J., Hook L, Pain R., Palanque- 
Delabrouille N., Perrett K., Regnault N., Rich J., Taillet 
R., Baumont S., Bronder J., Ellis R. S., Filiol M., Lusset 
V., Perlmutter S., Ripoche P., Tao C., 2006, ApJ, 648, 
868 


Suzuki N., Rubin D., Lidman G., Aldering G., Amanullah 
R., Barbary K., Barrientos L. F., Botyanszki J., Brodwin 

M. , Gonnolly N., Dawson K. S., Dey A., Doi M., Don¬ 
ahue M., Deustua S., Eisenhardt P., Ellingson E., Faccioli 
L., Fadeyev V., Fakhouri H. K., Fruchter A. S., Gilbank 
D. G., Gladders M. D., Goldhaber G., Gonzalez A. H., 
Goobar A., Gude A., Hattori T., Hoekstra H., Hsiao E., 
Huang X., Ihara Y., Jee M. J., Johnston D., Kashikawa 

N. , Koester B., Konishi K., Kowalski M., Linder E. V., Lu- 
bin L., Melbourne J., Meyers J., Morokuma T., Munshi 
F., Mullis C., Oda T., Panagia N., Perlmutter S., Post¬ 
man M., Pritchard T., Rhodes J., Ripoche P., Rosati P., 
Schlegel D. J., Spadafora A., Stanford S. A., Stanishev 
V., Stern D., Strovink M., Takanashi N., Tokita K., Wag¬ 
ner M., Wang L., Yasuda N., Yee H. K. G., Supernova 
Cosmology Project T., 2012, ApJ, 746, 85 

van den Bergh S., Li W., Filippenko A. V., 2005, PASP, 
117, 773 

Wang X., Wang L., Filippenko A. V., Zhang T., Zhao X., 
2013, Science, 340, 170 

White R. E., Keel W. C., 1992, Nature, 359, 129 

White HI R. E., Keel W. C., Conselice C. J., 2000, ApJ, 
542, 761 

Wood-Vasey W. M., Miknaitis G., Stubbs G. W., Jha S., 
Riess A. G., Garnavich P. M., Kirshner R. P., Aguilera 
C., Becker A. C., Blackman J. W., Blondin S., Challis P., 
Clocchiatti A., Conley A., Covarrubias R., Davis T. M., 
Filippenko A. V., Foley R. J., Garg A., Hicken M., Krisci- 
unas K., Leibundgnt B., Li W., Matheson T., Miceli A., 
Narayan G., Pignata G., Prieto J. L., Rest A., Salvo M. E., 
Schmidt B. P., Smith R. G., Sollerman J., Spyromilio J., 
Tonry J. L., Suntzeff N. B., Zenteno A., 2007, ApJ, 666, 
694 

Zibetti S., Chariot S., Rix H., 2009, MNRAS, 1410 


© 2010 RAS, MNRAS 000, 1-?? 



